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Abstract 

We calculate the 0{as) corrections to longitudinal spin-spin correlations in e^e^ — > 
qq. For top quark pair production the 0(as) corrections to the longitudinal spin- 
spin asymmetry amount to less than 1% in the g'^-range from above tt-treshold up 
to = 1000 GeV. In the e"^e~ bb case the 0{as) corrections reduce the 

asymmetry value from its m = value of —1 to approximately —0.96 for g-^-values 
around the Z-peak. This reduction can be traced to finite anomalous contributions 
from residual mass effects which survive the m — > limit. We discuss the role of 
the anomalous contributions and the pattern of how they contribute to spin-flip and 
no- flip terms. 
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Recently there has been renewed interest in the role of quark mass effects in the production 
of quarks and gluons in e~^e~ annihilations. Jet definition schemes, event shape variables 
and heavy flavour momentum correlations are affected by the presence of quark masses for 
charm and botton quarks even when they are produced at the scale of the Z-mass [0, Q . 
A careful investigation of quark mass effects in e'^e~ annihilations may even lead to a 
alternative determination of the quark mass values |l|, ^ 3 • There is obvious interest 
in quark mass effects for ti production where quark mass effects cannot be neglected in 
the envisaged range of energies to be covered by the Next Linear Collider (NLC). In the 
calculation of radiative corrections to quark polarization variables residual mass effects 
change the naive m = pattern of polarization results that results from the absence of 
spin-flip contributions in the m — > case [§, H, ^ H- 

In this report we calculate the 0{as) radiative corrections to longitudinal spin-spin 
correlations of massive quark pairs produced in e~^e~ annihilations. The longitudinal po- 
larization of massive quarks affects the shape of the energy spectrum of their secondary 
decay leptons. Thus longitudinal spin-spin correlation effects in pair produced quarks and 
antiquarks will lead to correlation effects of the energy spectra of their secondary decay 
leptons and antileptons. As a byproduct of our calculation we discuss the m —>■ limit 
and the role of the 0{as) residual mass effects. We delineate how residual mass effects 
contribute to the various spin- flip and no- flip terms in the m — limit. 

Let us begin with by defining the differential joint quark-antiquark density matrix dur = 
daxiX2;\[x'2 where Ai and A2 denote the helicities of the quark and antiquark, respectively. In 
this paper our main interest is in the longitudinal polarization of the quark and antiquark, 
and in particular, in their longitudinal spin-spin correlations. Thus we specify to the 
diagonal case Ai = A'^ and A2 = Ag. 

The diagonal part of the differential joint density matrix can be represented in terms of 
its components along the products of the unit matrix and the z-components of the Pauli 
matrix nx^ [nx-^ = Pi^ for the quark and (IT3 = for the antiquark, pi = Pi/\pi\)- One has 

rfjT = 1 (^dal O 11 + da^^^Vg ® ]1 + da^^^H ® trg + rfa^^^^^ V3 ® 4x3^ . (1) 

An alternative but equivalent representation of the longitudinal spin contributions can be 
written down in terms of the longitudinal spin components s{ = 2Ai and = 2X2 with 

s{, S2 = ±1 (or s{, S2 e {T> i})- One has 

da{s{, 4) = \{da^ da^'-\s[ + da^'^h'^ + da^'^'^\s[si) . (2) 

From CP invariance one knows that da and da^^'^^^^ obtain contributions from the 
parity-even VV- and y4A-current products, whereas da^^^^ and da^^^^ are contributed to 
by the parity-odd VA and y4\^-current products. The parity-even terms da and dcr'^^i^^) 
are C-even and thus symmetric under q ^ q exchange, whereas the parity-odd terms 
are C-odd and thus one has da^^^\pi,p2) = —da^^^\p2,Pi) for the single-spin dependent 
contributions. 
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Eq. (|^) is easily inverted. One has 



da = da 1) + da iU) + da il^) + da ill) 

da^'^^ = da{^^) + daiU) - da{l^) - daiU) 

da^''^ = dain) - daiU) + dail^) - daiil) (3) 

^ dain)-dain)-daii^) + daiU). 

0{as) radiative corrections to the rate component da have been discussed before ^ 
including beam polarization effects 0] and beam-event correlation effects [§, 0. The 0{as) 
radiative corrections to the longitudinal spin component da'^^^^ have been recently cal- 
culated including also beam polarization and beam-event correlation effects p. As 
concerns the longitudinal spin-spin correlation component da^^'^^^^ the 0{as) tree graph 



contributions have been determined in |T0[. Here we calculate the 0{as) radiative correc- 
tions to the fully integrated spin-spin correlation component cr(^i^2) -^^here we average out 
beam-event correlation effects. 

As before we write the electro-weak cross section e+e^ <l{Pi)<l{P2) and e+e~ 
(l{Pi)(l{P2)g{P3) in modular form in terms of two building blocks 0. Thus we write (beam 
polarization effects not included and beam-event correlations averaged out) 
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da{s{, S2) = ^ i^guidai + rfcrf '^'■'siSa) + gi2{da2 + da'i^'^^'' s[s2) 

+g,,{daf'h[ + daf'\s'S). (4) 

The index i = 1, 2, 4 on the rate components is explained later on. 

The first building block Qij = 1,2,4) specifies the electro-weak model dependence 
of the e~^e~ cross section. For the present discussion we need the components gu, gi2, gu, 
5(41, 5(42 and 5(44. They are given by 



gu = Qf -'^QfVeVfReXz + {Ve + (^e)i'^f + (^f)\Xz\ , 

gi2 = Q)-2QfVeVfReXz + {vl + al){v)-a))\xz\^, 

gu = 2QfVeafReXz -'2{vl + al)vfaf\xz?, 

g4i = 2QfaeVfReXz - '^Vetteivj + aj)\xz\^, (5) 

5(42 = 2QfaeVfReXz -^v^tteivf - a'j)\xz\'^ 

gM = -2QfaeafReXz + ^VetteVfaflx, 
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where, in the Standard Model, Xz{q^) = gM^q'^/{q^ - + iMzTz)~\ with Mz and 
Tz the mass and width of the Z° and g = GF(8v/27ra)~i ^ 4.49 ■ 10"^ GeV~^. Qf are 
the charges of the final state quarks to which the electro-weak currents directly couple; 
Vg and Oe, Vf and af are the electro-weak vector and axial vector coupling constants. For 
example, in the Weinberg-Salam model, one has Ve = —I + 4sin^ 9w, o,e = —1 for leptons, 
V/ = 1 — I sin^ 9w, a/ = 1 for up-type quarks {Qf = |), and fj = — 1 -|- | sin^ 9w, c^/ = ~1 
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for down-type quarks {Qf = — I). In this paper we use Standard Model couphngs with 
sin^ 9w = 0.226. 

The second building block is determined by the hadron dynamics, i.e. by the current- 
induced production of a quark- ant iquark pair which, in the 0{as) case, is followed by gluon 
emission. In the 0{as) case one also has to add the one loop-contribution. We shall work 
in terms of unpolarized and polarized hadron tensor components Hjj^l, -^i/+L' ^u+l ^^'^ 
Hjj^i where the spin decomposition is defined in complete analogy to Eq. (^. In the two 
body case e~^e~ qq the unpolarized rate components are given by 

- -^Hl.. (6) 

In the three-body case e^e~ qqg the unpolarized differential rate components and the 
unpolarized hadron tensor components Hljjj^^ are related by 

J^ = J^2Hh+Liy,^) (^ = 1,2). (7) 

As kinematic variables we use the two energy-type variables y = 1 — 2piq/q'^ and z = 
1 — 2p2<i/<f'- The same relations hold for polarized production. 

The index i = 1, 2 in Eqs. (P) and (|^) specifies the current composition in terms of the 
two parity-even products of the vector and the axial vector currents according to (dropping 
all further indices on the hadron tensor) 

1 1 

ttI _ (uVV j_ o-AAn tt2 _ ^ t ttVV ljAA\ (o\ 

In the parity-odd case one has 

1 

5' 

The notation closely follows the one in Thus the nomenclature {U + L) in Eqs. (^) 
and (|^ denotes the total rate {U: unpolarized transverse, L: longitudinal) after averaging 
over the relative beam-event orientation. 

The generalization of the above cross section expressions to the case where one starts 
with longitudinally polarized beams is straightforward and amounts to the replacement 

gu ^ [il-h-h+)gu + ih- -h+)gu] (10) 
gu ^ [{l-h-h+)gu + ih- -h+)g^,] (2 = 1,2) (11) 

where h~ and (—1 < h"^ < +1) denote the longitudinal polarization of the electron and 
the positron beam. 

Let us begin with by listing the Born term contributions to the various polarized and 



Ht. = 7;iHl^ + H^^)- (9) 



unpolarized two-body hadron tensor components. One has = 4m^/g^, v = a/1 — ^) 
Hlj^^iBorn) = (4 - ^Ncq', Hl^^{Born) = 3^Ncq', 



H]}it\BoTn) = -(4 - 30Ncq\ H'Jif\Born) = -iNcq\ (12) 

Hlj%{Born) = AvNcq\ H^%{Born) = -AvNcq^ 
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The 0{as) spin dependent hadronic tree-body tensor 

gluon spin 



(13) 



can easily be calculated from the relevant Feynman diagrams. The {U + L)-component 
is then obtained by contraction with the four-transverse metric tensor {—g^u + q^lv/'f)- 
Finally, the longitudinal spin components of the quark and antiquark can be projected out 
with the help of the respective longitudinal spin vectors. They read 



^(\/(i-i/)^-e;0,o,i-y) 



■'2 



{sir = ^(V(l-^)'-e;(l-^)sin^i2,0,(l-z)cos^^i2) 



(14) 



with 



cos 9 



yz + y + z-l + E, 
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(15) 



The resulting spin-independent and single-spin dependent components of the hadron tensor 
have been given before (see e.g. f^, H, ^). Here we list the spin-spin dependent piece. One 



has {vy :- 



'1 -y)^ :-- 







i(^ife) 

U+L 



VyVz 



4{i2-m+e) + {i-oi^-m 



1 1 

y2 ^2 



+ (4-30(8-70 (^ + ^]+2{12-50{y+z) 



y y 



z z 



-2(4-0(1/^ + ^^) -(4-30^^-^ + --- 

yz-^ z-^ y^ j 

-2(1 - 0(2 - 0(4 - 30- - (4 - 0(6 - 50 f- + - 

yz \z y, 

+2(4-50 (7 + 7) +4^^-^ 



(16) 



0-2(^1^2) 



VyVz 



-4e + (l-0e 



1 1 



z^ 



(8 - 70 



-6{y + z)- 2{y' + + 5 



z^ 



y Zj 



-2(2-0(1-0 



1 

yz 



(6+0 

\z y \z y^ 



-Ayz 



(17) 



What remains to be done is to perform the phase space integrations and to add in 
the one-loop contributions. In this calculation we perform the requisite two-fold phase 
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space integration over the full {y, z) phase space. As in ^ the infrared singularities 
are regularized by introducing a gluon mass. The infrared singularities in the tree graph 
and one-loop contributions cancel and one remains with finite remainders. For the sake 
of completeness we include in our results also the unpolarized hadron tensor components 
which are needed for the normalization of the longitudinal spin-spin asymmetry. The 0{as) 
corrections (tree plus loop) read 



1 

4' 



(4 - 0(2 - 0^ + 7(192 - io4e - + 



2(4 - (2 - 0(^8 - + 2viho + 2ti2; 



^(18-0^ + ^(24-8^-^^3 



-6 (2-0(t8-t9) + 2t;(tio + 2ti2; 



(19) 



0-1(^1^2)/ , 



N 



78C - 5^^ + 3^^) 40 + 32 - 14^ + 120/^ + + 3^%/^ 

2v ^ ^ 



16 - 42^ + 3ie- + 8(4 - 30v' ) % 



+2(4-30 (2 -0(^8 -tie) + 2t;(t 
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2t 



12) 



(4-0(8-3e-a^ 



-2(8 - lOe + e^)ti4 + (32 - 88^ + 76^^ - 19^^)% 

/1 .0 



(20) 



0-2(^1^2)/ , 



- (54 - 65^ + 3^^)— + 58 - 56^^ - 3^ - 3^^^ 

(J 

h 



+ 2 + e - 40 - 8t;^ K + 2 (2 - 0(^8 - tie) + 2t;(tio + 2t 



12, 



-(96 - 140^ + 350 -3C 



f211 



' 2(10 + 30ti4 + (8 - 20e + 130)% 



where we have used an overall normalization factor = asNcCpq^ /4:7rv. The unpolarized 
hadron tensor components Hlj_^_i^{as) and Hljj^]^{as) including the 0{as) rate functions U 
{i = 3, 8, 9, 10, 12) have been calculated before in 0. They are listed here for completeness. 
In addition to the rate functions calculated in |^ the spin-spin contributions bring in a set 
of new rate functions tj {i = 13, 14, 15, 16). The complete set of rate functions needed in 
the present application is given by 



In 



1+v 
1-v 
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^8 
^9 



l+V 

1-v 



= 21n 



In 



^10 
ti4 



= In 



+3 [Li2 



2v 



In 



+2Li2 



In 



1-v 
ti2 = In 



- Li2 

1 + z; 
1-v 

- Li2 



2f 



1 + 1; 

+ 2 ( Li2 
2v 



— Lio 



2v 



1-v 



l + v 



1 -V 



+ 



1 + t; 



ii3 = In 



1 + ^ ^ 



- 2Li2 + Li: 



1 - v 



— Lig 



l + v 



= In 



l + v 
1-v 



+ In 



^16 



In 



-4Li2 



ll-v 
l + v 

In 



.2-Vf, 

/ r 

+ 2Li2 ' ' 
Av^ ' 



^(1 + ^;)^ 



l + v 



+ 2Li2 

2v ^ 



2-v?r 



1. 



+ -Li2 
2 



:i+v)^ 



Li2 

2 



(1 + vyj 

{l-v)\ 



-2v 



{l-vf 



(22) 



We are now in the position to discuss the normahzed longitudinal spin-spin correlation 
function (P^^) which is defined as 



^ ^(TT) + ^(n) + ^aT) + ^(u)' 



(23) 



The mean (P^^) is taken with regard to all phase-space variables including the beam-event 

orientation variables. 

Before we present our numerical results we briefly pause to discuss the m — > limit (or 
equivalently the — > oo limit) of the various hadron tensor components and the longitu- 
dinal spin-spin correlation functions. In the m — > limit the hadron tensor considerably 
simplifies and one has (-f^^+i — > 0) 



H-i (J J\ - \(lJ^ I 0-1(^1^2) / A 



= Ncq'[{l-s{4)[l+^-^C, 



4 3as 



Cfs[s2 



3 47r 



(24) 



By hindsight the 0(0;^) contribution has been split into a normal piece and an anomalous 
piece (indicated by square brackets). The normal piece proportional to (1 — sjsa) is defined 
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spin 
config. 


vv 




m = 


m ^ 


m = 


m ^ 


(TT) 





t = o+[t] 








(Ti) 


2 




2 


l = 2-[|] 


(it) 


2 




-2 


-| = -2 + [|] 


(ii) 





l = o+[t] 









Table 1: ©(a^) corrections to specific spin configurations in QCD(m = 0) and QCD(m 
0). Entries are given in terms of contributions to the hadron tensor components 
Hi;^i{si,S2) and S2) in units of SagCpNcq'^ /4:Tc. Anomalous contributions are 

shown in square brackets 



in terms of the usual massless QCD(m = 0) no-fiip 0{as) result whereas the anomalous 
piece is a residual mass effect which survives the m — limit. 

Let us complete our discussion of the m — > limit by also listing the corresponding 
m ^ limit of the longitudinal single-spin hadron tensor components. One has 



j(^^c/+L'Si + iifz+LSgJ {2b) 

s{-si). (26) 



2 3as„' 

— ■ O F 

3 47r 



In Table |l| we list all m — contributions to the various spin configurations where again 
we have split off the m = no-fiip contributions. 

At threshold the production cross section is dominated by the s-wave vector current 
contribution. One thus has Hlj_^^i^{s{, s^) = H1jj^j^{s{, $2) and Hljj^j^{]]) = |if^^j;^(ti) = 
^if^^_^(J.t) = i7fy_|_^(||). From Eq. (p3|) one then obtains a threshold value of (P^^) = 
-1/3. 

Let us now present our numerical results. In Fig. 1 we plot the mean longitudinal spin- 
spin asymmetry (P^^) against the cm. energy ^/q^ for top quark pair production. The 
longitudinal spin-spin asymmetry rises from its threshold value of (P^^) = —1/3 to around 
(P^^) = —0.9 at y/q^ = 1000 GeV. The 0{as) correction to the asymmetry amounts to 
less than 1% in the g^-range from above ti treshold to \/q^ = 1000 GeV. 

In Fig. 2 we present our results on (P^^) for bottom quark pair production starting from 
bb threshold (where {P") = -1/3) up to \/q^ = 100 GeV. For the lower g^-values from 
threshold to about 30 GeV the 0(0;^) corrections are quite small. Starting at around \/q^ = 
30 GeV the 0{as) correction become larger. The Born term contribution very quickly 
acquires its asymptotic limiting value (P^^) = — 1 due to the fact that the corrections 
to the leading term are quadratic in the ratio m/y/q^. Contrary to this the 0{as) curve 
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remains below the naive limiting value of —1. From the limiting formula 



1+^ 

7r 




V 37r J 


U 


a 


s 



TT 



one concludes that a large part of the deviation is made up by the anomalous contributions. 
For example, at the position of the Z pole the limiting value of the anomalous contribution 
to (P^^) amounts to {P^^){ano'm) = 0.048 {as{mz) = 0.118). From the full calculation one 
finds {P^^){Born) = -0.996 and = -0.964. Thus the deviation of (P^^) from 

its naive value of (P^^) = — 1 can be seen to arise to a large part from the anomalous 
contribution. 



Note added in proof: When preparing this manuscript for publication we became aware 



of a preprint on the same subject by M.M. Tung, J. Bernabeu and J. Penarrocha [pT 
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Figure Captions 

Fig. 1: Energy dependence of 0(1) and 0{as) mean longitudinal spin-spin correla- 
tions (P^^) in e~^e~ — > ti{g) 

Fig. 2: Energy dependence of 0(1) and 0{as) mean longitudinal spin-spin correla- 
tions (P^^) in e~^e~ bb{g). The vertical line indicates the 66-threshold 
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